Positive current cross-correlations in a highly transparent normal-superconducting 
beam splitter due to synchronized Andreev and inverse Andreev reflections 
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Predictions are established for linear differential current-current cross-correlations dS a ,b/dV in 
a symmetrically biased three-terminal normal metal-superconductor-normal metal (NSN) device. 
Highly transparent contacts turn out to be especially interesting because they feature positive 
dS a ,b/dV. At high transparency, processes based on Crossed Andreev Reflection (CAR) contribute 
only negligibly to the current and to dS a ,b/dV. Under these circumstances, current-current cross- 
correlations can be plausibly interpreted as a coherent coupling between the two NS interfaces in 
the form of synchronized Andreev and inverse Andreev reflections (AR-AR), corresponding to the 
process where a pair of electron-like quasi-particles and a pair of hole- like quasi-particles arrive from 
the normal electrodes and annihilate in the superconductor. Hence, positive dS a ,b/dV does not 
automatically imply CAR. For tunnel contacts, dS a ,b/dV is positive because of CAR. In between 
these two extremities, at intermediate transparencies, dS a ,b/dV is negative because both processes 
which cause positive correlations, occur only with small amplitude. We use scattering theory to 
obtain analytic expressions for current and noise, and microscopic calculation using a tight binding 
model in order to obtain a clear interpretation of the physical processes. 



I. INTRODUCTION 

Beautiful experiments on transport and noise in nor- 
mal metal-superconducting (NS) hybrids allow probing 
the microscopic physics associated to the superconduct- 
ing condensate and quasi-particles. For instance, the 
superconducting gap A was revealed by tunnel spec- 
troscopy on a normal mctal-insulator-superconductor 
tunnel junction. In NS structures with highly transpar- 
ent contacts, Andreev reflection 1 is the phenomenon by 
which pairs of electron-like quasi-particles from the nor- 
mal electrode N can enter the superconductor S and join 
the condensate. Additional processes appear in N a SNt 
structures with two normal electrodes N a and N&: an 
electron coming from N Q may be transmitted as an elec- 
tron into Nft (elastic cotunneling, EC), or it may be 
transmitted as a hole into Nj, (crossed Andreev reflec- 
tion, CAR 2-36 ). The amplitudes of these two processes 
decrease exponentially with a characteristic length scale: 
the coherence length £, which is inverse proportional to 
the energy gap A in the ballistic limit. Therefore, three- 
terminal nanoscale devices with distance R > £ between 
the contacts are especially interesting: conductance and 
noise experiments on them probe both the condensate 
and the quasi-particle states. 

Concerning current-current correlations, theorists have 
envisioned two kinds of experiments for the long term: 
using entanglement in quantum information devices, and 
testing entanglement in the electronic Einstein- Podolsky- 
Rosen (EPR) experiment 4 ' 25 ' 37 The EPR experiment is 
not considered here, but instead the basic problem of 
current-current cross-correlations 38 in NSN structures is 
addressed. Some experiments based on NISIN structures 
have been reported recently' Our task here is not to 
understand the tunnel limit, where an insulating oxide 
layer I is inserted in between the normal and the super- 
conducting electrodes, but the opposite limit of highly 



transparent interfaces where Coulomb interactions^ are 
not expected to play a predominant role. 

Based on the limiting case of tunnel contacts^' 
one may erroneously conclude that positive differential 
current-current cross-correlations dS a ,b/dV are equiva- 
lent to CAR. However, this is not the case because, as 
we show, positive dS a ^/dV can well be obtained in the 
absence of CAR. 

Current-current cross-correlations are negative 4 ^' for 
non-interacting fermions. A flux of bosons leads to pos- 
itive cross-correlations* 3 - and negative cross-correlations 
are found for bosons impinging one by one onto a beam 
splitter. 44 Cross-correlations can be positive in interact- 
ing fermionic systems' - ' as well as in multiterminal 
NS structures" 

The recent experiment' and other experiments under 
way involve normal electrodes separately connected to a 
superconductor, 28 with a geometry similar to that con- 
sidered in the following. Not only the noise can be eval- 
uated in various set-ups, but also the full histogram of 
the charge transmitted in a given time interva h 25 ' 27 i 30 
In addition to these set-ups, relevant information is 
also obtained from "zero dimensional" chaotic cavities 
in contact with a superconductor, in connection with 
the possibility to observe positive current-current cross- 
correlations due to CAR' 

In what follows, attractive interaction binding pairs 
of electron-like quasi-particles is present everywhere in 
the superconducting region. Two electron-like quasi- 
particles of a pair injected into the superconducting re- 
gion remain glued by the BCS mean field interaction, 
in contrast with the dissociation of a Cooper pair enter- 
ing a chaotic cavity (see Ref. [23h ■ The physics behind 
current-current cross-correlations in the NSN structure 
with highly transparent contacts considered here was not 
really elucidated in Ref. l28l . in spite of the important ob- 
servation that dS a .b/dV is positive but the non-local con- 
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ductance is negative at low bias. Unusual properties can 
be realized with the following experimental conditions: 
First it is assumed that the same voltage V a = Vb = V 
is applied on the two normal electrodes N a and Sec- 
ond the temperature is very low, and third high val- 
ues of interface transparencies are used. Direct elec- 
tron transmission is Pauli blocked at zero temperature 
because of the same voltage applied on the two normal 
electrodes. On the other hand, crossed Andreev reflec- 
tion is strongly reduced at high transmission, regardless 
of the applied voltages V a and Vb- These two features 
of EC and CAR are shown to be a sufficient condition 
for deducing dS a ,b/dV > 0, not due to CAR, but due 
to what is called here synchronized Andreev reflection 
and inverse Andreev reflection (AR-AR). The positive 
current-current correlations at high transparency are not 
in disagreement with what is found in Ref. for strong 
coupling to the superconductor, where a gap is induced 
in a chaotic cavity. The cross-over between high val- 
ues of interface transparency and tunnel contacts will 
also be investigated: CAR has a dominant contribution 
to dS a ,b/dV for tunnel contacts, and dS a fi/dV is nega- 
tive at the cross-over for intermediate values of interface 
transparency, where AR-AR is suppressed. 

The article is divided into two independent main sec- 
tions where current-current correlations are evaluated on 
the basis of (i) the scattering approach for a homogeneous 
superconducting gap (see Sec. and of (ii) microscopic 
calculations taking into account the strong inverse prox- 
imity effect (see Sec. IIVI) . Technically, these two sections 
rely on different approaches and are based on different 
assumptions. Overall agreement between the two calcu- 
lations is obtained. Final remarks are provided in Sec. El 
Technical details are as much as possible left for Appen- 
dices. The article starts with a preliminary section con- 
taining definitions and a summary of the main results. 



II. PRELIMINARIES 
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FIG. 1: (Color online). The two models studied in this paper: 
Panel (a) shows the one-dimensional geometry used in the 
BTK calculation, which is a simplified description for a NSN 
structure in three dimensions. The central superconducting 
electrode has length R and the superconducting gap Ao is 
uniform. Panel (b) shows a representation of a three-terminal 
device described by a tight-binding Hamiltonian. The two su- 
perconducting reservoirs have the same phase and thus con- 
stitute a single terminal. 



where a is the projection of the spin on the quantization 
axis and the sum over n runs over the M tight-binding 
sites describing the interface. Tight-binding sites on the 
normal side of the interface N a S are labeled by a n , and 
their counterparts in the superconducting electrode are 
labeled by a n . The hopping amplitudes between elec- 
trode N a and the superconductor are denoted by t an an 



A. Current, noise and current-current correlations 

We start with general definitions of current and 
current-current cross-correlations. The geometry of the 
considered set-up is shown in Fig. [TJ The central island 
is superconducting, and it is connected by highly trans- 
parent contacts to the two superconducting reservoirs on 
top and bottom. The superconductor S is made of the 
central island and of the two reservoirs. 

The operator giving the current flowing at time t from 
the normal electrode N a to the superconducting island S 
at the N a S interface [see Fig. [TJb)] takes the form 

M 

/„, a (t)= (^n,a„C+„,a(i)c a „, CT (i) (1) 

<T,n— 1 

(t)Ca»,*(t)) - 



of a magnetic field. The average current I a = (I a ,a{t)) 
is the expectation value of the current operator given by 
Eq. ©. 

Current-current auto-correlations in electrode N a are 
given by 

SaAf) = (si a (t + t r )si a (t)) + (si a (t)6i a (t + t')), (2) 

with 5I a (t) = I a (t) — (I a (t)). In the absence of ac- 
excitations, the average current given by Eq. ([1]) is time- 
independent and the auto-correlations S a , a (t') given by 
Eq. ([2]) depend only on the difference t' of the time ar- 
guments. 

Similarly, current-current cross-correlations between 
the electrodes N a and N& are given by 

s a ,b(t') - (si a (t + t')si b (t)) + (si b (t)6i a (t + t')), (3) 

where h describes the current at the interface with the 
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FIG. 2: (Color online). Two independent local AR processes 
are shown on panel (a), for two interfaces separated by a dis- 
tance R much larger than the coherence length £. The labels 
1 and 2 on the diagram correspond to electrons and holes re- 
spectively. For each AR process, an electron impinging from 
the normal electrode onto the interface is converted into a 
hole, and a pair is transmitted into the superconductor. If 
R ~ £, the two AR processes can be coupled coherently by 
non-local propagation in the superconductor [see panel (b)], 
with a quartet as an intermediate state and penetration of 
a charge 4e into the superconductor, which is qualitatively 
equivalent to double CAR. Panel (c) shows independent AR 
and AR processes at the two interfaces, supposed to be far 
apart. Panel (d) shows AR-AR for R ~ £, which is qualita- 
tively equivalent to double EC. 



normal electrode N&. Zero- frequency auto-correlations 
(S a ,a) and cross-correlations (S a ,b) are defined as the 
integral over t' of <S a ,a(*0 and of S a , b (t'). The defini- 
tion used here for the Fano factor is as follows: F a a = 
S a , a /2el a and F ab = S a , b /2ey/I a I b , with I a = I b . With 
this definition, the Schottky formula leads to F = q*/e 
for a Poisson process transmitting a charge q*. 



B. Known results and open questions 

A few facts related to the non-local conductance for 
arbitrary values of interface transparency are already 
known. Only local Andreev reflection AR and local in- 
verse Andreev reflection AR come into account for con- 
tacts separated by a distance much larger than the co- 
herence length [see Figs.^a) and (c)]. Local AR means 
that an electron is converted into a hole and a pair of 
electron-like quasi-particles is transmitted into the super- 
conductor, and local AR means that a hole is converted 
into an electron and a pair of hole-like quasi-particles is 
transmitted into the superconductor, which annihilates a 
Cooper pair. Local AR and AR contribute to transport 
if the separation between the interfaces is much larger 
than £, and in two-terminal configurations where the su- 
perconductor is not connected to ground. 

Other quantum processes appear in a three-terminal 
configuration if the distance between the contacts is com- 
parable to the coherence length: CAR and EC. However, 
non-standard types of "non-local" processes can be also 
obtained by merging AR at the interface N Q S to AR at 
the interface SN&, forming what is called here AR-AR 
[see Fig. C^b)]. Physically, AR-AR would correspond 
to the synchronized transmission of two pairs from elec- 
trodes N a and N b into the superconductor, which can be 
seen as double CAR. However, this process does not con- 
tribute to non-local transport at zero temperature. Con- 
versely, AR at interface N a S might be associated to AR 
at interface SN b , leading to AR-AR. The corresponding 
non-local resistance is independent on the value of inter- 
face transparency in the tunnel limit i 12 ' 13 Qualitatively, 
AR-AR can also be seen as double EC. 

The non-standard non-local process AR-AR in- 
volving pairs appears naturally when expanding 
diagrammaticallyi^ the non-local conductance to order 
i 8 , with t the hopping amplitude at the interfaces. As 
it is shown below, AR-AR plays a central role in under- 
standing the positive^ current-current cross-correlations 
for highly transparent contacts, in a regime which is not 
described by perturbation theory in t. 

A very recent preprint^! points out the possibility 
of "synchronized Andreev transmission" in the current- 
voltage characteristics of a SNS junction array. Synchro- 
nization manifests itself in this work as specific features 
in the current-voltage characteristics of the two-terminal 
SNS junction array. We arrive here at the conclusion 
that synchronization of Andreev processes is also possi- 
ble if the separation between the NS interfaces is com- 
parable to the coherence length. The dominant chan- 
nel AR-AR is shown to result in positive current-current 
cross-correlations. 

As mentioned in the introduction, experiments on 
current-current cross-correlations in NSN structures have 
already started^ A few basic questions regarding 
current-current cross-correlations for highly transparent 
contacts have not yet received a satisfactory explanation. 

First what is the physics behind the positive^ 8 - linear 
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differential current-current correlations dS a .b/dV for a 
highly transparent NSN beam splitter? It is shown that 
dS a ,b/dV > is not an evidence for CAR (which would 
prevail 26 for tunnel contacts in the absence of Coulomb 
interactions). An interpretation in terms of AR-AR is 
proposed. 

Second, how do current-current cross-correlations de- 
pend on the sample geometry? Current-current cross- 
correlations decay with the geometry-dependent coher- 
ence length, as it will be obtained from microscopic cal- 
culations in Sec. IIVI 

Third, what is the value of cross-correlations at in- 
termediate transparencies? Experimentalists can realize 
tunnel or highly transparent contacts, by oxidizing or not 
the sample during fabrication. Intermediate values of in- 
terface transparency are more difficult to control but it is 
nevertheless useful to quantify how "perfectly transpar- 
ent" the NS contacts should be in order to obtain AR-AR. 
A cross-over from positive to negative dS a ,b/dV is found 
in BTK calculations as the normal state transmission co- 
efficient Tjy is reduced below a value typically of order 
T N ~ 1/2. 

Fourth, do the predictions established with a ballistic 
superconductor hold also for a disordered superconduc- 
tor? It is shown at the end of Sec. HVI that. for strong in- 
verse proximity effect, AR-AR is responsible for positive 
cross-correlations also in the case of a disordered super- 
conductor in the regime where the elastic mean free path 
is shorter than the coherence length. 



III. HOMOGENEOUS SUPERCONDUCTING 
GAP: BTK CALCULATION 



Latin labels fc, / run over a, b, referring to the two nor- 
mal electrodes N a and N^. Greek labels a,f3,-f,8 denote 
electrons or holes in the superconductor. The notation 
/ e (cj) = 9(eV — Li) stands for the distribution function of 
electrons at zero temperature, and fhiy) = 0(—eV — u>) 
is the one of holes, where 6(x) is the Heaviside step- 
function. In Eqs. (|4]) and <j5j) , sgn(a) = +1 if a = e, 
and sgn(a) = — 1 if a = h. 

The elements of the s-matrix are evaluated from the 
BTK approach (see Appendix [A} for a one-dimensional 
N a SNb junction [see Fig. QJa)]. Step-function varia- 
tion of the superconducting gap at the interfaces is 
assumed. A repulsive scattering potential V(x) = 
H [S(x) + 5(x — R)] is introduced at the interfaces. The 
transparency of the interfaces is related to the BTK pa- 
rameter Z = H/hvp, with the Fermi velocity Vp. The 
interface transparency is characterized by the value of the 
normal state transmission coefficient Tjv = 1/(1 + Z 2 ). 
Highly transparent contacts correspond to Tn = 1, and 
tunnel contacts correspond to Tjv <C 1. 

In the one-dimensional model considered in this sec- 
tion, current and noise are highly sensitive to the length 
R of the superconducting region: they oscillate as a 
function of R with a period equal to the Fermi wave- 
length Xp <C R. These oscillations can be interpreted 
as Friedel oscillations where the contacts with the nor- 
mal electrodes play the role of impurities. In a more 
realistic higher-dimensional model with more than one 
transmission mode, the oscillations in the different modes 
are independent and thus are averaged out. In order to 
simulate qualitatively multi-dimensional behavior with a 
one-dimensional system, current and noise are averaged 
over one oscillation period&S^ 
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The BTK approach^ allows calculations of the 
current- voltage characteristics of a NS point contact with 
arbitrary interfacial scattering potential. It was first gen- 
eralized in Ref. [ID and later in Ref. H3 to the case of 
non-local transport. Useful physical informations can be 
obtained even though the gap is not self-consistent in the 
BTK calculation. We will consider a one-dimensional ge- 
ometry, as shown in Fig. [If a). 

The current li at the interface with the normal 
electrode Nj, and the current-current correlations S% t j 

between the two normal electrodes N, and N, are can be neglected and sff(u) « sff(u = 0). Using this 



Sfj(R) = T-l drS id (r). 

AF JR-Xf/2 



(7) 
(8) 



In the studied limit of small applied voltage eV -C A, 
the energy-dependence of the scattering matrix elements 



expressed^ in terms of the scattering matrix s as 



h 



J ' dui sgn(a)sgn(/3) 



i>M (4) 

(5) 



approximation, it is possible to perform these integrals 
analytically. 

In this limit, the current I{ through the normal lead 
Nj is obtained as 



2y/2e 2 VT 2 (2 - T N ) cosh(i?/g) sinh(i?/g) 2 
h [2(2 - T N ) 2 cosh{R/£) 2 - 8(1 - T N )]* 



(9) 



with 



The local current-current correlations Su at the inter- 
face NjS and the current-current cross-correlations Sjj 
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FIG. 3: (Color online). Variations of the linear differential conductance [left panels, in units of e 2 /h], linear differential auto- 
correlations [center panels, in units of e s /h] and linear differential cross-correlations [right panels, in units of e 3 /h] as a function 
of the normal interface transparency Tjv, for the values of R/i indicated in the figures. A single-channel one-dimensional BTK 
calculation is used [see Fig. [Ha)]. The lower Panels show that dl% v /dV and dS^ a /dV are almost independent on R/i for 
R/i > 1, and that dS^ v b /dV becomes very small as R/i increases above ~ 1. The insert of the bottom right panel shows the 
variations of the normalized current-current cross-correlations [dS^ b /dV(TN)]/[dS a ,b a ^ /dV(TN = 1)]. The data corresponding 
to R/i = 2, 3, 4, 5 superimpose after this rescaling. 
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FIG. 4: (Color online). Variations of the Fano factors F a , a = S^u/lel^" [top panels] and F a ,t = S%j,/2el% v [bottom panels] as a 
function of the normal interface transparency Tjv, for the values of R/i indicated in the figures. A single-channel one-dimensional 
BTK calculation is used [see Fig. [Ha)]. The ratio between the gap Ao and the Fermi energy ef is Ao/ef = 10~ 4 in this 
simulation. The insert of bottom right panel shows the variations of the normalized crossed Fano factor [F a ,b(TN)]/[F a ,b(TN = 
1)]. The data corresponding to R/i — 2, 3, 4, 5 superimpose after rescaling. 



between the interface N^S and the interface SNj are 
6AV2e 3 VT^ smh(R/e) 2 sech{R/e) 

i i 7 

h(2 - TV) 3 [2(2 - TV) 2 cosh(i?/£) 2 - 8(1 - TV)] 5 
- 128(1 - TV) 4 + 2(2 - TV) 2 (1 - TV)(96 - TV 

x (192 - TV(116 - TV (20 - 37V)))) cosh(i?/£) 2 
- (2 - T N ) 4 (72 - TV(144 - 7V(82 - T N 
x (W + T N ))))cosh(R/£) 4 

+ 8(2-T Ar ) 6 (l-r Ar )cosh(i?/e) 6 

(10) 



MV2^VT% sinh(fl/£) 2 sech(fl/£) 
h{2 - TV) 3 [2(2 - TV) 2 cosh(i?/£) 2 - 8(1 - TV)] 5 
128(1 - Tat) 4 - 2(2 - T N f{l - T N ){2 + TV) 

x (2 - 3TV)(8 - (8 - T N )T N ) cosh(i?/£) 2 

+ (2 - TV) 4 (8 - TV(16 + T N (2 - TV(10 + T N )))) 

x cosh(R/0 4 ■ 

(11) 
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The linear conductance dl£ v /dV, the auto-correlations 
dS%" a /dV and the cross-correlations dS%" b /dV given by 
these formulas are shown in Fig. |3] as a function of the 
normal transmission coefficient T^, for different values 
of the distance between the contacts. As it is expected, 
the conductance increases with interface transparency for 
R/£ ^ 1- The conductance depends on the ratio 
while R is smaller than the coherence length £, but it does 
almost not vary as R is increased above £. For R/£ < 1, 
the conductance is non-monotonous when plotted as a 
function of (see Fig. [3] top left). An explanation for 
the non-monotonous behavior is the enhanced transmis- 
sion due to the finite size of the superconductor^ 

Starting from tunnel contacts, the linear differential 
auto-correlation dS^ v a /dV first increases with interface 
transparency as larger current leads to larger noise. The 
differential noise reaches a maximum and almost vanishes 
for perfect transparency if R > £, as it is expected for 
a single NS junction. 38 Differential current-current cross- 
correlations dS® v b /dV are positive for R/£ < 1 in the 
extreme cases of very high and very low interface trans- 
parency, and take negative values in between, as shown in 
Figs, fright), and in the insert of Fig. [^bottom right). 

The Fano factors F a>a = S*™/2e/™ and F afi = 
S lv b /2el av correspond to the noise normalized to the cur- 
rent, which allows to get rid of the trivial effect that, at 
low transparency, the noise increases when the current 
increases. The variations of F a<a and F a ^ feature local 
minima at intermediate values of T/y. In the insert of 
Fig. @|bottom right), the Fano factor F a ^ is shown for 
different R/£ > 1. As it can be seen, for R/£ > 1 the 
Fano factor is independent of R/£ after normalizing to its 
value for T/v = 1. The Fano factor F a>b is positive outside 
the region of the minima. For R/£ > 1, F aM takes the 
value F a a ~ 2 for T/v -C 1, and the value F a a ~ for 
T N = 1. ' 

We first make some remarks in order to confirm the va- 
lidity of our calculation. Expected behavior is recovered 
in some known limiting cases: 

(i) For tunnel contacts (Tjy C 1), the Fano factors are 
given by F a , a « 2 - ±sech(i?/£) and F a ,h sa ±sech(i?/£), 
which is in agreement with the expected limiting values 
F a , a = 2 and F a>b = for R/£ — > oo. This corresponds to 
the doubling of the effective charge for Andreev reflection 
at a single NS interface in the tunnel limit i 52 1 53 

(ii) For highly transparent interfaces (T/v ps 1), the 
Fano factors vanish exponentially with increasing as 
F a . a — F a ,b « sech(i?/£), as it is expected for a single 
highly transparent NS interface. 

(iii) dS% v b /dV and F ayb are positive and very small for 
R/£ > 1 and T/v — 1, in agreement with Ref. [26. Only 
CAR contributes to current-current cross-correlations for 
T N < 1. 

(iv) As it can be seen in Fig.^top left), the linear con- 
ductance dl%° jdV is suppressed for R/£ < 1 because the 
number of normal states within the gap energy is ~ R/£ 
for this geometry. For a three-dimensional grain it is pro- 
portional to ~ kpR 3 /£, which can be much larger than 



unity even for R/£ < 1. The suppression of Andreev 
processes at high transparency for R/£ <C 1 is thus not 
expected to occur in the case of a three-dimensional su- 
perconducting grain. 

The Schottky limit is realized for low values of interface 
transparencies, which lead to ({SI a — 8I b ) 2 ) av — Ael%" '. 
On the other hand ((SI a + 8I b ) 2 ) av = 4e(I™ + Ig v ). One 
concludes that F a>a ~ 3/2 and F ayb ~ 1/2, in agreement 
with the plateau obtained in Figs. HJtop) for the depen- 
dence on T N of the Fano factor. Eqs. ®, <jT0]) and (fTTj) 
reproduce these values in the limit R/£ — > and T/v <C 1. 

For highly transmitting interfaces TV = 1, we ob- 
tain F a ^ a = F a ,b- These values can be confirmed by a 
more simple calculation in the limit R/£ -C 1 (see Ap- 
pendixE]). The identity F a , a = F a>b implies that I^-i^ 1 ' 
is noiseless for T/v = 1, independent of i?/£: 



dT< («/„(<) - 5I b (tj) (8I a {t + r) - SI b (t + r))) ( 



/ qciv I qav n qciv \ 

~ 2 V°a,a ~r O b ^ b — £O a b ) 

jav 

= ~y (F a ,a + F b , b - 2F a>b ) = 0. 

By comparison, in a fermionic beam splitter with highly 
transparent contacts, a charge e transmitted from the 
source to N a means no charge transmitted from the 
source to N b . Thus, for fermions, it is the sum I a + I b 
that is noiseless. 

We can see from Fig.[3]that the linear differential cross- 
correlations of current noise dS%" b /dV are positive for 
T^v -C 1 and Tjy — 1, while they take negative values 
in between these two limiting cases. For T/v -C 1, the 
positive dS% v b /dV can be explained by the presence of 
CAR processes. This is in agreement with perturbative 
calculations carried out in Ref. l26l However, for perfectly 
transmitting interfaces T/y — 1, CAR processes do not 
occur, as the elements of the scattering matrix describing 
CAR (e.g. s^ b ) equal zero for T/v = 1 (see Appendix IB) . 
In this Appendix we obtain in the limit R/£ <C 1 



Joan 
g3 a,b 

dV 



and 



3CAR = 



,e,ft 



for R/Z < 1 
= T EC for R/£ > 1 



(13) 



= 0. 



(14) 



In section IIV1 we will use a microscopic model in order 
to obtain an understanding of the processes contributing 
to the noise, and to explain dS" v b /dV > in the absence 
of CAR. 

To summarize this section, we have obtained analyt- 
ical expressions for the current and noise in the one- 
dimensional BTK model. It was shown that dS% v b /dV > 
0, and that Eq. (JB) holds for R/£ > 1 and for T N ~ 1. 
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Cooper pair splitting dominates for small normal trans- 
mission coefficient <C 1, while what will be inter- 
preted in Sec.|lV]as AR-AR dominates for Xjv — 1. These 
two processes are suppressed for intermediate Xjy, result- 
ing in negative cross-correlations in this parameter range. 



IV. MICROSCOPIC CALCULATIONS 

The current-current cross-correlations for highly trans- 
parent interfaces can be further investigated in the two- 
dimensional tight-binding set-up shown in Fig. [TJb). 
First, using analytic calculations, we analyze the differ- 
ent contributions to the noise S a ,b and show the absence 
of contributions due to CAR processes. The positive con- 
tributions are attributed to processes which we refer to 
as AR-AR (this notation stands for Andreev reflection 
and inverse Andreev reflection) , a denomination which is 
motivated by the microscopic analysis. Second, numeri- 
cal calculations are performed, which take into account 
the inverse proximity effect by determining the gap in 
a self-consistent manner. In addition, disorder will be 
included in the calculation. As in the last section, we 
restrict our study to the case where the same voltages 
V a = Vb = V are applied on both normal leads and V is 
small compared to the gap A. 

The expression of current-current cross-correlations 
S a ,b can be decomposed as a sum of six contributions 
according to the types of transmission modes in the su- 
perconductor: 



S a 



Si 



CAR 



s 



EC 



s 



AR-AR 



Sar- 



AR 



s' 



Smixed, 

where the expression and the meaning of the different 
contributions are provided in Appendix |TJ] Scar con- 
tains the noise attributed to crossed Andreev reflections, 
which contains transmission modes in the electron-hole 
channels. Sec, t ne noise due to elastic cotunneling, con- 
tains transmission modes in the electron-electron or hole- 
hole channels. With AR-AR, we refer to synchronous lo- 
cal Andreev reflections at both interfaces, while AR-AR 
links a local Andreev process at one interface to a local 
inverse Andreev process at the other one (see Fig. |2J). 

The Andreev-reflection is highly localil (compare with 
the BTK calculations in Sec. IIII1 where it vanishes ex- 
actly). This motivates the assumption |G a '^| <C \G^\, 
which is confirmed by our numerical calculations also in 
the presence of the inverse proximity effect (see Fig. [5J . 
Using this simplification, one obtains Scar = Sar-AR = 
0. The terms S' and Smixed vanish after integration 
over the energy uj. Thus the total current-current cross- 
correlation S a .b depends only on the term S^r-AR- 

In order to understand what type of microscopic pro- 
cesses are described in AR-AR, we start from the formula 
giving the current-current cross-correlations in terms of 



0.003 



11.1)1)2 



0.001 




FIG. 

.-,1,2 

Gi' 1 



d = \tln ~ Vm\ 

5: (Color online). The anomalous Green function 
Jm is small compared to the normal Green function 

3 . As the elements G an ,p m (between site n on the left 



interface and site m on the right interface) decay exponen- 
tially with d = \yn — TJm\, we have evaluated the mean values 
\(Gct n ,0 m \}n,m of the elements with the same value of d. The 
main contribution to non-local transport comes from small 
values of d. While the main frame shows the results obtained 
for a constant gap A, the inset shows the corresponding data 
with the self-consistent gap. 



the Keldysh Green functions G + ' and G ,+ : 



S a ,b{w) =Tr 



6,0 u a,p 



G+ 5 - - + 

P,a U a,b 



(16) 



where the trace is carried out over the Nambu labels and 
the different transmission modes at the interfaces. It can 
be shown that all terms in S^r-ar are obtained from 
the anomalous contributions of the type G +,_ ' ' 2 G~' + ' 2 ' 1 
and G +, ~' 2,1 G~ ,+ ' 1,2 . Let us consider one of these terms 
as an example (the same conclusions are obtained for all 
terms) and suppose again that G a '^' 2 = if a and fi are 
on different interfaces (that is, crossed Andreev reflection 
does not contribute), and that a symmetric bias voltage 
is applied on the two normal electrodes. One has 



G 



+,-,2,1, 



b,a (t,t')G-^(t',t) 
'- (Ca4(Oc 6 ,t(i))(c+ f (r 



a, J, 



(*))■ 



(17) 



This equation can be understood as a relation between 
initial and final states, as it is shown in Fig. [6] In gen- 
eral, these states can be connected by many different pro- 
cesses. However, the microscopic formula for the current- 
current cross-correlations [see Eq. (ID4[) ] shows that the 
initial and final states are related by an Andreev process 
at interface N Q S, and, at the same time, by an inverse 
Andreev process at interface SN& [see Fig. [Bfb)]. In an 
Andreev process, an electron is converted into a hole and 
a pair of electron-like quasi-particles is transmitted into 
the superconductor. In an inverse Andreev process, a 
hole is converted into an electron and a pair of hole-like 
quasi-particles is transmitted into the superconductor. 
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FIG. 6: (Color online). Schematic representation of how 
AR-AR couples to current-current cross-correlations. The ini- 
tial state consists of (i) a spin-down electron created at a at 
time t, on the superconducting side of the N a S interface, and 
of (ii) a spin-up electron destroyed at time t at b, on the 
normal side of interface SNt. The final state consists of a 
spin-down electron destroyed at time t' at a, on the normal 
side of interface N a S, and of a spin-up electron created at 
time t' at /3, on the superconducting side of interface SNb. To 
the AR-AR process shown on Fig. 01, the permutation of two 
fcrmions is added. Taking into account the resulting minus 
sign leads to positive current-current cross-correlations. 



The pair of electron-like quasi-particles annihilates with 
the pair of hole-like quasi-particles and the remaining 
electron and the hole are exchanged between the two in- 
terfaces. This results in dS AR _j^/dV > 0. 

In addition to the analytic calculation, we performed 
numerical simulations in order to analyze a more realistic 
model. Details about the used method^- are presented 
in Appendix[C] The self-consistent simulations presented 
below take into account the inverse proximity effect cor- 
responding to the reduction of the superconducting gap 
within a distance ~ £ from the contacts. Self-consistency 
is equivalent to current conservation for the electrons in- 
jected from the normal reservoirs and transmitted into 
the superconducting ones. 

The values for the linear differential cross-correlations 
dS a ,b/dV , and for the EC and CAR transmission coef- 
ficients Tec and Tcar are plotted as functions of the 
length N of the superconductor in Figs. [7] and |U The 
notations Tec and Tcar refer to the transmission modes 
in the superconductor (advanced-advanced or retarded- 
retarded modes not exchanging electrons and holes for 
Tec , and advanced-retarded Green functions exchanging 
electrons and holes for Tcar)- Tec and Tcar are given 
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FIG. 7: (Color online). The figure shows dS a ,b/dV as 
a function of N for highly transmitting interfaces [for N 
see Fig. D^b)]. In agreement with the BTK calculations 
dS a ,b/dV > 0. The exponential decay is described by the 
coherence length which increases with M, in agreement with 
Ref. H3. The values M — 20 (red squares) and M — 25 
(blue circle) are used. Strong deviations from the expo- 
nential decay (red lines) appear for N ~ M (see Ref. I20T ). 
The data points are obtained from Eq. (|C4p . The decay is 
ocexp[-2ATao/£(M)]. 
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FIG. 8: (Color online). The exponential decay of the EC 
transmission coefficient Tec [Eq. (|18[) . open red squares for 
M = 20 and open blue circles for M — 25] , and of the crossed 
Andreev transmission coefficient Tcar [Eq. (|19[) , filled green 
squares for M = 20 and filled purple circles for M = 25] is 
shown. The decay is oc exp[-2Na /£(M)], with £(20) = 8ao 
and £(25) = 9a . 



by 



Tec 
Tqar 



W 2 Tt 



( c a,i,i 



G 



R,l,l 



cc,/3 



R,l,2 



(18) 
(19) 



where W is the hopping amplitude in the bulk and at 
the interfaces, and a n and f3 m run respectively over all 
the sites on the superconducting side of the N a S and 



SN), interfaces. The notation G 



A(R), 



stands for the 



Nambu component (m,rij) of the advanced (retarded) 
Green function connecting i to j. 

The dependence on M of the coherence length was 
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already found in a previous work<2H The differential cross- 
correlations dSa^/dV are positive, which is in agreement 
with the preceding BTK calculation. Differential cross- 
correlations dSa.b/dV show exponential decay (see Fig. [7]) 
as a function of N, because the two normal electrodes 
N a and N& are coherently coupled by evanescent states 
in the superconductor. The BCS coherence length as 
obtained from the fits fulfills R x /£ > 1 in a wide range 
of simulation parameters. This is the range in which the 
BTK calculation leads to dS^JdV = 8^T EC . 

For a highly transparent NS contact, one has 



G 



2,1 



G 



/3m, £>„ 



1/4 W 2 where W is the hopping am- 
plitude in the normal and superconducting electrodes, 
and at the interface (see Appendix [E]). The identity 

T ar A ar - T ^ holds if l G ^l « 1^1 (see Ap- 
pendixlEj, with 



r A,A 

" AR-AR 



-W 2 Tr {G^Gf'f 



°a,/3 U /3,a 



(20) 



where the superscript "A, A" refers to an advanced- 
advanced transmission mode where electrons and hole 
are conserved. With these assumptions, the total noise 
can be written as 



£ = 3.8ao, 3.6ao for V/W — 1.00, 1.25, as compared to 
£ = 8eto with V/W = in the ballistic limit. The coher- 
ence length in the presence of disorder becomes smaller 
than its ballistic value, as for a superconductor in the 
dirty limit. Within error-bars, Eq. ([21]) is fulfilled also 
in the presence of weak disorder, while the coefficient in 
Eq. (|2"2")l is changed, resulting in 



,°' b — 40— Xcar- 
dv h 



(24) 



It is concluded that Tcar <C Tec implies that S a ,b — 
S^-^ and that Eq. j21j is fulfilled. Thus, it was shown 
that, in this parameter regime, dS a ,b/dV > is evidence 
for AR-AR, not for CAR. 

The BTK approach and the microscopic calculations 
lead to positive dS a ,b/dV > 0, which is due to the ex- 
change of fcrmionic quasi-particles (see Fig. [HI . This is 
the main physical result of our article: dS a ,b/dV > at 
high transparency is not due to CAR. It is due to pairs 
of electron-like quasi-particles, pairs of hole-like quasi- 
particles, and exchange of fermions. 



dV 



i—T A ' A 

h AR-AR 



rriA.A rr 

T ar-ar- Tec - 



(21) 
(22) 



CONCLUSIONS 



Considering the numerical data, plots of dS ai b/dV as 
functions of Tec and ?car (see Fig. [9]) , and a comparison 
between T A V A and T EC confirm Eqs. {HJ and (|32J and 

A xx- A xx 

show in addition 



-W * AR ' 



(23) 



Comparing with the previous BTK calculation, it is 
suggested that Eq. (f23|) is model-dependent, because 
Tcar = for the BTK model while Tcar is finite but 
small in the self-consistent microscopic calculation. On 
the other hand, Eq. (l2"Tj) , which is obtained also for the 
BTK model, is expected to be a fundamental relation. 

Until now, we only studied perfect systems without dis- 
order. The numerical calculations based on microscopic 
Green functions give us however the additional possibil- 
ity to add disorder to the model, which is an unavoidable 
ingredient to describe experiments. 

Eq. (|21[) remains approximately valid in the presence 
of weak disorder (see Fig. I10[) . introduced in the form 
of a random on-site potential uniformly distributed in 
the interval [— V, V], with elastic mean free path l e ~ 
(W/V) 2 ao smaller than the ballistic coherence length ~ 
(W/A)ao, but not small compared to the Fermi wave- 
length. The parameter W is the hopping amplitude in 
the bulk of the superconductor, taking the same value at 
the interfaces because of highly transparent contacts. 

The coherence length is reduced as the strength of 
disorder increases. The coherence lengths are fitted to 



The article was already summarized in Sec. Ill Bl and 
thus we conclude with a brief overview and final remarks. 
We have evaluated current-current cross-correlations in 
a NSN structure with a homogeneous superconductor 
(without self-consistency in the order parameter), and 
with strong inverse proximity effect (with self-consistent 
microscopic calculations for a two-dimensional three- 
terminal set-up). For both approaches, the linear differ- 
ential cross-correlations dS a ,b/dV are positive for highly 
transparent contacts and decay exponentially with a 
characteristic length set by the coherence length. Pos- 
itive dS a ,b/dV arises in this set-up not because of CAR, 
but because of what is identified as the correlated pene- 
tration of pairs of electron-like quasi-particles and pairs 
of hole- like quasi-particles into the superconductor in the 
form of AR-AR. The positive sign of dS a ,b/dV is due to 
the additional exchange of two fermions. It is emphasized 
that the proposed mechanism does not involve quartets in 
the superconductor, because AR-AR does not contribute 
to the current-current cross-correlations. Direct evalu- 
ation of dS a ,b/dV leads to dS a ,b/dV = 4(e 3 /h)T^ A -—, 

A xl- A Jrv. 

and to — ~ Tec, which holds also for a super- 
conductor with weak disorder and with elastic mean free 
path shorter than the coherence length. 

Finally, correlations between pairs of Andreev pairs 
were discusse d 59 ! 60 in connection with noise in an An- 
dreev interferometer. 

Future studies of this class of set-ups might include 
the effect of arbitrary applied voltages and quantitative 
effects of three-dimensional geometry. 
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Tec Tcar 

FIG. 9: (Color online). Relation between dS a ,b/dV and the transmission coefficients Tec and Tcar: Panel (a) shows the 
relation between dS a ,b/dV and Tec and panel (b) shows the relation between dS a ,b/dV and Tcar- The data are the same as 
in Fig. [7] The different points correspond to different lengths N along Ox axis [see N in Fig.QJb)]. The blue dashed lines show 
the linear fits dS a ,b/dV = 8T EC , and dS a ,b/dV = 50T C ar. 
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FIG. 10: (Color online). Relation between dS a ^/dV and the transmission coefficients Tec and Tcar in the presence of weak 
disorder (V/W = 1 and V/W = 1.25): Panel (a) shows the relation between dS a ,b/dV and Tec and panel (b) shows the 
relation between dS a ,b/dV and Tcar- The figure is obtained by averaging over the contributions between the combinations 
of all different transverse sites x n and x m in the two normal leads. The width of the model is given by M = 20 [see M in 
Fig. Eft)]. 
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Appendix A: Details on the scattering approach 

The elements sf'- of the scattering matrix are calcu- 
lated within the BTK approach^ using two-component 
wave-functions describing electrons and holes respec- 
tively. The indices i,j refer to the normal electrodes N a 



and Nfc, while a, f3 run over the two components describ- 
ing the electrons and holes. 

For example, the wave-functions for an electron incom- 
ing from electrode N a take the form 

rf> a (x)= Q (e^ + sZe-^*) (Al) 



4>s(x) = 



ipb{x) = 




(A3) 

where ip a (x), ips(%) and ipb(x) are the parts of the wave- 
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functions in the electrodes N OJ S, and N& respectively [see 
Fig. Ufa)]. The notations q^ + \ q(~', k^ + ' and stand 
for the wave-vectors in the normal and superconducting 
electrodes: 



fc (±) =k F ± »/£ 
q {±) =k F , 



(A4) 
(A5) 



h.e 



and s b ' e c can be deter- 



with > 1. 

The elements s e a f a , s£;~, s b 
mined using the continuity of the wave-functions at the 
interfaces [tp a (0) = ^s(O) and ips(R) = ">Pb(R)] an d the 
boundary condition for the derivatives [*p' s (0) — ?p' a (0) = 
Hi> a (0) and ip' b (R) - i>' s (R) = H^ b {R)]. The BTK pa- 
rameter Z is defined by Z = H/Kvf- 

The remaining elements of the scattering matrix can 
be obtained from the other possible scattering processes 
(e.g., a hole incoming from electrode N&) by analogous 
calculation. 

By comparing the different equations, the symmetry 



(A6) 




of the scattering matrix can be obtained. 

The BCS coherence factors uo, vq appearing in these 
equations are given by: 



(A7) 



The coherence factors uq and vq are interchanged under 
complex conjugation and changing sign of the real part of 
energy (a small imaginary part is supposed to be added 
to uj). 



Appendix B: Current-current cross-correlations for 
the BTK model in the limit Z ~ 



For highly transparent interfaces (Z ~ 0), the expres- 
sions obtained for sf'f simplify considerably, local re- 
flections (s^f = 0) and non-local Andreev reflections 



(sf'j — 0, with i ^ j,a ^ /3) are suppressed^ Thus, 
only local Andreev reflection (sff, with a ^ 0) and 
transmission without branch-crossing (s"'", with i =/= j) 
can occur. The non-zero elements of the scattering ma- 
trix are given by 



-R/i 



-R/i 



(Bl) 
(B2) 



Assuming s{uj) to be constant in the range [—eV,eV] 
(we study the case eV <C A), Eqs. (UJ) and J3]) can be 



written as 

4e 2 

T - — V s e ' ft 
4e 3 



(B3) 



S a , b = -^\V\ (s^stptlXy + s h a fjsl%stp h b: ?) 

(B4) 

for i ^ j. In the limit R/£, <C 1, the elements of the 
scattering matrix (evaluated for uj — > 0) are given by 



S b,b ~ 1 



R 



h,h 



leading to 



> 0. 



(B5) 
(B6) 

(B7) 



The result for S a , a is obtained by an analogous calcula- 
tion. 

Combining Eq. (fB3| for I a with Eq. (|B7)) for S a . a and 
Sb,b leads to F a , a ~ F b ,b ^ 1 for eV <C A , Tn — 1 and 



Appendix C: Technical details on microscopic 
calculations 



The microscopic calculations are based on the following 
tight-binding Hamiltonian on a square lattice: 



Rs = ~W ( C n.<7 ( + <7„ V 



n 

^ 1 ^n c n, cr c n,cn 



(Cl) 



with the hopping amplitude W between nearest neighbor 
sites n and m separated by a distance ao- The normal 
electrodes are described by an analogous Hamiltonian 
with no pairing term and no disorder. Highly transparent 
contacts with interfacial hopping equal to W are used in 
Sec. lIVI The parameter A„ is the superconducting order 
parameter at site n. It is determined self-consistently in 
Sec. lIVI on the basis of the recursive algorithm developed 
in Ref. I2fj. The gap in the superconducting reservoirs 
takes the fixed value Ao. Disorder is introduced at the 
end of Sec. IIVI in the form of a random on-site poten- 
tial V n on each tight-binding site, uniformly distributed 
in the interval [— V, V]. The gap and phase profiles in 
the superconducting island are shown in Fig. Qj] in order 
to illustrate the output of the part of the code perform- 
ing the self-consistent calculation. Because of disorder, 
the gap fluctuates strongly from one tight-binding site 
to the next. The phase profile shows a smooth exponen- 
tial decay from the NS interfaces. The accuracy of the 
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(a) Gap profile (b) Phase profile (semi-log scale) 




FIG. 11: (Color online). The figure shows the gap profile in panel (a), the phase profile in semi-log scale in panel (b), for a 
given realization of disorder with V/W = 1. See the text for the meaning of the parameters V and W. The x and y axis are 
the same as in Fig. [T^b). This figure is obtained with the self-consistent algorithm developed in Ref. [53 for the ballistic case. 
The phase in panel (b) is small but non-zero because the algorithm takes into account non-equilibrium effects. 



self-consistent calculation gives access to variations of the 
phase over almost two orders of magnitude. 

The average current in Eq. (TT]) is evaluated from the 
Keldysh Green function^ 



In 



2e 



(C2) 



Zero frequency noise (see Sec. [TTJ) is obtained 
from S a .a = 2e 2 t 2 /h J dujS a . a (uj) and S a . b = 
2eH 2 /h fdu)S a , b (u)), with 



S a ,a («)=£>[ G+;- an {uj)G- a + am («) 
n ' m ~+Gi- Q >)G-,;+>) 

-G+- a >)G-+,>) 
-G+;; 0n (a;)G-+ m (a;)" 



(C3) 



(CM) 



The trace is evaluated over the Nambu labels. Eq. (|C4j) 
is a generalization of Ref. 56 to two interfaces with many 
channels. The numerical calculations presented in the 
main body of the article are based on Eqs. t|C2|) . (IC3|) 
and (ICil). 
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Appendix D: Microscopic calculations for the noise formula 

In this Appendix we provide the complete formula for the noise S a ,b(co) given by Eq. (|C4j) for sub-gap voltage 
< eV < A). This allows considerable simplification of the expression of S a .b because the Keldysh Green function 
g H = in the isolated superconductor, as there exist no single-electron states. The total noise can be decomposed 
into different terms according to the types of transmission modes in the superconductor as: 

S a ,b(w) — Scar + Sec + ^ar-ar + Sar-ar + Smixed + S'. (Dl) 

In the normal lead N, the Keldysh Green functions read g^/ 11 = 2i-Kp^np(uj — eV), g^ly 22 — 2iTrp^np(uj + eV), 
Sattv' 11 = 2iitppfnp(—L0 + eV) and g~^ 22 = 2iTrpNnp(—ui — eV), where p^ is the density of states at the interface in 
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the normal lead, and, at zero temperature, tif(x) — 6(—x), with 6(x) being the Heaviside step-function. This gives 

_i ^ ]_ \- 22 A 22 h 11 

9n n 9n n = 9n n 9n n = 0> lading to further simplification of the expression for S ,b. 

The contribution Scar is given by the advanced-advanced or retarded-retarded transmission modes in the electron- 

a 12 2 1 

hole channel, in the form of the combinations of the type G Q ^ G S . Physically, these microscopic processes can be 
interpreted as Cooper pair splitting as appearing in CAR. The expression for Scar i s as follows: 



— hH„+- ,22 / ^H,2l / ^fl,l2/ 1 , nir <A,22.2 „ \/i i njMUJ. \ 
bb 9aa G af3 G p a I 1 + 2lG ac nt Pa ) (1 + 2lG^ ntp b ) 

r,-nA,22 ,2 \ 
2lG^ TTt p b ) 



Scar = 2£ 4 1^ + g 66 y aa 

. -+,22+-Al r ,R,12 r <R,2l/ 
+ 9bb 9aa <J Q( 3 ^ p a \ 



9aa 



> pa v ' eta r< 

,22 +-4ir< A >2i.^A,i2/ ■ or<-R,22 .2„ \/ ; , orWi ( 2„ \ 

G a,3 G ,3a + 2G aa p a )(l + 2& p ^Pb) 

f a 2 \i + 2G^ 1 nt 2 Pa )(i + 2G«£\t 2 p b ) 



9aa 
-+.!!„+ 



,22^,12^,21 /. , O( ofi,l] 



(D2) 



where t is the hopping amplitude at the interfaces. 

The contribution Sec contains advanced- advanced and retarded-retarded transmission modes in the electron- 

_A_ 12 _A 21 

electron or hole-hole channel, in the form of the combinations of the type G a ^ G a' a . This is the contribution 
to the noise of normal electron transmission in the form of EC. The expression reads 



Sec = 2t 



+ 9bb ' 9. 



G::rG': : rii-2G^^ Pa ){i-2G^\ep b 

2G^ 2 irt 2 p a )(i-2G^ 22 irt 2 p b 



9bb 



— h,22 



,22 r »fl,22 r ,ii,22/. 
U ap U Pa \ l 



g^g^G^G^Hi + 2G^nt 2 p a )( l + 2G^nt 2 Pb ) 



? ii,ii 



(D3) 



' 22 G A J 2 G A p f^ + 2G£" TrfpaXi + 2G^nt 2 Pb ) 



iR,22 



«,22 ,2, 



Now we consider processes that do not appear in lowest order perturbation theory. First the contribution S^r-ar 
contains advanced-advanced and retarded-retarded transmission modes, in the form of the combinations of the type 
Ga'^Gpa" 2 ■ Physically, this can be interpreted as the contribution to the noise of synchronized Andreev and inverse 
Andreev processes (AR-AR), with the exchange of two fermions at the same time, thus leading to dS a: b/dV > (see 
Fig- EJ. This contribution dominates the current-current cross-correlations dS a ,b/dV in the considered set-up. This 
contribution reads 



C _ 0^2,8 „ / , _ — h,22 + -,n r A.12 r yA,21 r ,R,22 r ,RSl - + .11 „+- ,22 n A,21 n A,12 n R,ll n R,22 

AR-A"R - _S7r 1 Po-Pb I + 9 bb 9aa ^aa ^ /3B '-'ap U Pa + 9 b b 9aa ^ aa ^ pp ^ a p U pac 

(D4) 



-+,ll / ,+-,22 r tA,ll r YA,22 r ,fl,21 r <it,12 , n -+,22 n +-,U r iA,22 rl A,ll r iR,12 r iR,21 
a 9bb ^ap ^Pa U act ^ pp ' 9aa 9bb U a.p U 'pa U «a ^0P 



Another contribution not appearing in lowest order contains processes involving electron-hole conversion both at 
the interfaces and during propagation in the superconductor. The transmission modes in the superconductor are 
of the type G a 'p G ^ ' a . These terms correspond to the synchronization of two Andreev processes (AR-AR) . The 
corresponding expression reads 



C — 8^2,8 n ( , - + ,22 +-,22 r <A,21 r <A,21 r ,R,12 r ,R,12 - + ,11 ,11 n A,12 n A,12 n R,21 n R.,21 

"AR-AR — »7T t p a Pb I + 9 bb 9aa '-'aa (j pp a p °"/3c* + 9bb 9aa « aa " pp (j a p ^Pa 

, — h.22 H — .22 n A,2\ n A, 2\ n R,\2 n R, 12 — \-,U H — ,11 n A,12 ^A.12 n R,21 n R,21^ 
T 9aa 9bb U aP U Pa U aa ^pp +9aa 9bb U aP U pa U aa U pp 
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Another contribution contains transmission modes of the type G A,1 G A ' 12 . These terms read 



S' =8ir 2 t s p a p b 



i +-,11^,11 / -+,22^,21^,12^,11 
9aa ^ao [ilbb u /30 U af3 U pa 

, „+-Al n AA2 ^-+,22^,22^,22^,21 
+ 9aa ^oa \9fc6 ^,3/3 <J a 

-+-,22.^4,21 / - + ,ll r ,A,ll r ,fl,ll r ,fl,12 
9aa '-'aa [9bb U f3f3 °a/3 U pa 

+ -.22 r A,22 ( - + ,ll n A,12 r ,R,21 r ,R,22 
9aa u aa \9bb U f3f3 U af3 U f3a 



— h 11^,12^,^,11^^,2 
9bb U l3/3 ^a/3 U /3a 



') 



+ 



+ 9 a 



- + ,ll r ,A,ll r< R,21 l ^R,ll\ 

9bb U /3f3 U al3 U /3a ) 

- + .22 r ,A.22 r ,R,l2 n R,22\ 
9bb U pP U a/3 U l3a ) 

— K22 ^4.21 - 
366 



U PP U al3 U /3a 



G 



A,21 ( +-,11^,11^,12^,11 
a£ \9fc6 ^oa "/S/S 

i - + .22™4,22 / +-,11^,12^,22^,21 
+ 9aa ^a/3 ^66 °/3« b oa ^,3/3 

, -+,ii n A,ii („+-' 22 n A ' 21 r iR ' 11 r i 

+ 9aa ^a/3 \9bb U (3a ^ aa ^ 

4- n -+ A1 r A < 12 (n+-> 22 r A < 22 r R > 21 r R '' 

' 9aa '-'ap [9bb u (3a ^ aa ^0/3 
i- ( 6„ [ , „-+,22^fl,12 ( +-,11^,21^,11 

AtiTt p a [ + g aa G a ^ [g bb G aj3 G (ia 



+ -.22 r A.22 r R.22 n R^2 

9bb ^/3a ^aa ^pp 

+-,22^,21^,12^,22 
9 b b U pa U aa Lr, 



ft, 12 



22 



+-,11^,12^,21 
966 U ffa U aa 



) 

U PP J 

') 



966 



+-,11^,11^,11^,2 
9bb '-'pa '-'aa ^/3/3 



,22 n A,22 n A,21 
'-'ap ^fia 



9bb 



_1_ „- + .llr , -R>21 / +-,22^,12^,22 
+ 9aa '-'aa ^66 b o/3 

, _+-, 11/^4,12 / - + ,11^,21^,11 -+ 
+ 9aa ^aa ^9&6 U a/3 U Pa ~ 9bb 



,ll r< A, 11 r , A, 12 
22 / -,ft,22 ( ^ ( ft,21 N 



+ 4iirt b p b 



, „+-,22 r A,21 / - + ,22 r K,lz r i 
+ 9aa ^aa ^66 b a(S 

, -+,22,-yA,21 / +-,22^,22^,12 
+ 966 °/3/3 ^9aa '-'ap ^ f) a 

, -+,11™4,12 / +-,ll rY ft,ll r ,ft,21 
"T 966 ^0/9 [9aa '-'ap '-'pa 

. + -,22^ft,12 ( -+,11^,11^,21 
+ 9 6 6 ^^/S 1^9aa ^ aj 3 ^ f3a 

, +-,11^,21 / - + ,22^,22^,12 
+ 966 °/3/3 ^9aa U f3a 



9bb 



9aa 



9aa 



+ ,11qR,11qR,12 



Al r ,R,12 r ,R,ll 
°a/3 '-'Pa 



G 



~9a 



R,21 n R,22\ 
a/3 '-'fia J 



22 r iA,21 n A.22 
U aP U pa 

,11^,12^,11^ 



(D6) 



The last term involves advanced-retarded transmission modes 



S 



MIXED — 



87T 2 t 8 



'PaPb 



- AM 6 p a 



! H — .ll/ rl A.ll r ifl.l2 _ r iA,12/-Yii,22\/ / -y^,21 r< iJ,ll _ n A,22 n R,21\ 
9aa y^act u aa ^ aa u aa JV-* a(3 U pa ° ap U pa I 

ir<A.2iriRAi r<A,22r<R-2i\{r< A11 r< R12 r < A i X2 n R > 22 \ 

y^aa ^aa ~ U aa ° aa ll^aP U Pa ~ U a p U pa > 



- + ,22 n +-,ll(r<A,ll 

"9, 

-q-+- 11 q+-- 22 
■yaa iiaa 

— h.22 H — .11 
966 966 

— h,ll H — ,22 



+ 966 966 \ U pa U 



, n A,U n RA2 
\ U pa ^aP 

>A,21 n R,ll 

'ap 



r ,AA2 r ,R,22^ rr ,A,21^,RAl 
■ <J- Q( g ){^t 



ls-lt 

T pp u pp 



n A, 22 n R,21 \(siA, 11 
" ^o/? 



G 



RA2 



PP ^pP 



a/3 ^ 
,A, 22^,^,21 

^-,.4,12^^.22^ 
'-'PP U pP > 



-G 



— hll H — ,22/,^, 
966 966 ( G 

— h,22 H — .11 
— a,. a.. 



A, 12 
,A,21 



n R.12s (r <A, 1\ n RA\ 
~ U PP ){ U Pa U aP 



n A,22 n R,21-> 
U Pa U aP I 



,-2 H — ,11/^,21 . ^,21% ln A, l\ n RA2 n AA2 n R,22^ 
- 9bb 966 ( U PP + U PP )( U Pa U ap ~ U Pa U aP ) 

M^rfi n \ -l „-+, 22 „+- A1 (r iAA2 _L n RA2 \(n A ' 21 r iR > 11 r< A ^ 2 r< R ^ 21 \ 

-AlTTt P6[+9oa 9aa (^aa + G aa ){G a0 G 0a - G a0 G p a ) 

~~i n A.21 , ^,21^/^,11^,12 ( r-,.4,12 rr ii,22x 
\S-aa +^aa ){^ a p ^ p a ~ U ap U Pa > 



- a -+- 11 a + - a2 (^ A ' 21 

ifaa Vac 



(D7) 



Appendix E: Evaluation of the Green function G a ,/3 connecting two interfaces 



It will be shown that, at small energy compared to the gap, the advanced-advanced AR-AR transmission mode 
(^a'p' 1 ^'p'a 2 ) av can b e replaced by the opposite of the advanced-retarded EC transmission mode — (G^'^' 1 G^'^' 1 ) a «- 
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The Green function is expanded in powers of the exponential coefficient exp(— R x /£) appearing in g a ,B, giving 8 

G£p = MAjApNfrp + O {{gi P f) (El) 
Mi* = (j- 9iJa, a g£X*) ~* (E2) 

nf,p = {i- kbdtAA) ~* ■ (E3) 

This expansion leads to 



(A, 1,1 A, 2, 2 . ■ ( A,l,2 , A,2,l\ ■/ A, 1,1 , A,2,2\ , A, 1,2 A, 2,1 

9 a ,p -9 a ,p + l \9 a ,B +9 a ,p ) *(fl a ,> +g a ' tP )+9 a ',0 -9 a ] p 

■ ( A, 1,1 , A,2,2\ . A, 2.1 A, 1,2 A, 1,1 . A. 2, 2 , ■ / A, 1.2 . A, 2,1 

l {9 a ,3 + 9 a ,B )+9 a ,p -9 a ,p -9 a ,p + 9 a ,B + l {9 a ,p + 9 a ,B 

One has a = g^B = g a ,B for energies within the gap. The off-diagonal Nambu components are vanishingly small if 
w -C A. It is deduced that 

(G^Gf^' 1 }^ = {(g 1 ^ - g 2 a 2 p ) )av + {{g^ji + g 2 a ^) } av , (E6) 



leading to the identification of the advanced-advanced AR- AR transmission coefficient to the opposite of the advanced- 
retarded EC transmission coefficient, for small energy compared to the gap and for Nao > £. Eq. (IE5[) has no imaginary 

part because {9^ a9a b) av = at small energy compared to the gap. 

The fully dressed local Green function can be evaluated approximately by inverting the Dyson equation G QiQ = 

g a , a + g a ,ai a ,a9a,Ja,aG a , a , leading to G^'i' 2 = G^' 1 = 1/2W at energy small compared to the gap. 
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